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ABSTRACT 

We demonstrate ground-based submillimagnitude (< 10^"^) photometry of widely 
separated bright stars using snapshot CCD imaging. We routinely achieved this pho- 
tometric precision by (1) choosing nearby comparison stars of a similar magnitude and 
spectral type, (2) defocusing the telescope to allow high signal (> 10^ e~) to be ac- 
quired in a single integration, (3) pointing the telescope so that all stellar images fall 
on the same detector pixels, and (4) using a region of the CCD detector that is free of 
nonlinear or aberrant pixels. We describe semiautomated observations with the Super- 
nova Integrated Field Spectrograph (SNIPS) on the University of Hawaii 2.2m telescope 
on Mauna Kea, with which we achieved photometric precision as good as 5.2 x 10~^ 
(0.56 mmag) with a 5 minute cadence over a 2 hr interval. In one experiment, we mon- 
itored eight stars, each separated by several degrees, and achieved submillimagnitude 
precision with a cadence (per star) of ~17 minutes. Our snapshot technique is suitable 
for automated searches for planetary transits among multiple bright stars. 

Subject headings: Stars, Extrasolar Planets, Astronomical Techniques 



1. Introduction 

High-precision stellar ph otometry is used to dete ct planets orbiting ot her stars, exploiting 



the p henomena of a transit (ICharbonneau et al.l |2000| ) , a secondary eclipse (jCharbonneau et al 



20071 ) , or a microlensing event (jBond et al.l |2004| ) . Although photometry from s pace can be vastly 



more precise - Kepler achieves ~ 2 X 10"^ in 30 minutes with V ~ 10 stars (jKoch et al.l l20ld ) 
- the preponderance of astronomical resources are on the ground and must contend with absorp- 
tion, scattering, and scintillation by the atmosphere. Ground-based searches have been able to 
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detect Neptune- to Jupiter-sized planets on close-in, transit ing orbits around solar-mass stars (e.g., 
Henry et al.ll2000bl: iPollacco et al.l 120061 : iGillon et al.l 120071 ) , a multi-Earth-mass planet around an 
M dwarf (ICharbonneau et al.l 12003 ). and planets on more distant orbits around M stars: (e.g., 
Gould et ablboiol l. 



Relative photometric stability of a few millimagnitudes on the timescale of minutes to hours 
is sufficient to detect the transit of a Jupiter-sized planet (Fig. [1]). However, superior stability 
is required to resolve degeneracies between the unknown impact parameter, size of the star, and 
limb-darkening effects, and to search for transit timing variation indicative of other, unseen planets 
(jWinn et al.ll2009l : ISouthworthll2010l ). Furthermore, high-pr ecision Doppler studies are now detect- 
ing Earth- to Neptune-mass planets o n short-period orbits ([Mayor et al.ll2009l : iHoward et al.ll2009l : 
Lo Curto et al.ll2010l : IVogt et al.ll20ld ). The transit depth of these objects, even around smaller M 
stars, is at most a few millimagnitudes (Fig. [T]). Signal-to- noise ratio (S/N) can be improved by 
phasing data with the known or hypothetical orbital period (jHolman et al.ll2006l : IWinn et al.ll2007l : 
Gillon et al.l 120071 ). but the gains from phasing (and binning) are limited by time-correlated (red) 
noise, which has a (1//)'^ varyi ng power spectra l density, and does not decrease as the square root 
of the number of observations ( Pont et al.ll2006l ). Rather, photometric stability better than 10~^ 
(~1 mmag) in a single observation is required. 



To surpass a precision of 10 



TO^ e ^ must be acquired per integration. However, the 



dynamic range of optical CCD detectors limit the number of electrons in a single pixel to ~10^. 
One method to avoid exceeding the linear range is to read out the de tector at a rate of many hertz 
and to co-add the signal from multiple reads with low read noise (jSouza et al.ll2006l ). Another 
is to spread the signal over hundreds of pixels by sh a ping the point-spr e ad function (PSF) w ith 



orthogonal charge transfer ( Howell et al. 20031. 2005 : Tonry et al 



20051 : 



Johnson et al 



2009) 



or 



by defocusing the telescope (e.g., iGillon et al.l l2009l : IWinn et all l2009l : ISouthworth et all l2009l ). 

Longer integration times reduce scintillation noise and lower the fraction of time spent reading out 
the detector. 

The deleterious effects of the atmosphere can be reduced by performing differential photometry 
with one or more comparison stars observed in the same field of view as the target star. That correc- 
tion's accuracy will be limited by the signal from the comparison star (or collection of comparison 



stars). This technique commonly achieves photometric sta bility better t 



ran 1 mmag on ground 



based observations of variable stars and exoplanet transits ( 


Moutou et al. 


2004; 


Tonrv et al. 


2005: 


Hartman et al. 2005: Johnson et al. 2009: Gillon et al. 2009: 


Winn et al. 


2009; 


Southworth et al. 


2009, 


2010). 



A larger telescope can collect more signal but, all else being equal, its detector will have a 
smaller field of view with which to include a suitable comparison star. For a given /-number, the 
field of view of a telescope with a specified imaging detector is inversely proportional to the aperture 
D. The number-fiux relation for nearby stars has an approximate slope of —1.5, thus for a given 
integration time, the signal from the brightest comparison star in the field of any target star scales 
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as X or D^^^ (not considering binaries). In principle, there is a modest advantage with 

larger telescope aperture. However, for a typical CCD field of view of 10', the probability of an 
equally bright comparison star appearing in the same field becomes la rge (i.e., 80% at a galactic 
latitude of 30°) only by an apparent magnitude of 13 (|Zombeckl 120071 ) . This brightness is at the 
limit of current Doppler techniques, and would exclude those systems most amenable to follow-up 
observations. 

If a more suitable comparison star can be found in a separate, but nearby field, observations of 
the target and reference can be obtained by slewing the telescope between the two positions, and the 
exposure time can be adjusted to acquire equal signal from each star. The comparison star could 
also be the target of a transit search. These sequential observations will be neither simultaneous 
nor through the same column of the atmosphere. This technique has been employed for decades 
with single-channel photomultiplier tubes (PMTs) and has been applied to studies of (and searches 
for) e xoplanet trarisits and stellar variability, commonly taking advantage of automatic telescope s 



(e.g., I Genet et al 



1987 



Young et al 



1991 



Henrvl 119991 : iHerirv et al.M2000al : iGrauer et al.l l2008l V 

One of the seminal articles on the subject. lYoung et al. ( 199ll ). describes the limits of photometry 



when slewing between stars and making use of high-speed robotic telescopes. The authors suggest 
that the limit of g round- b ased autor nated photometry is ~1 mmag when all sources of error are 
properly handled. iHenrvl (jl999l ) and iHenrv et al.l (j2000al ) use automatic photoelectric telescopes, 
all less than 1 m in aperture, and achieve photometric precision >1 mmag for single observations 
by observing a set of approximately four stars at a time, of which three are comparison stars. They 
claim their precision is limited more by variability of comparison stars than by noise inherent to 
the technique. So far, none of these efforts have been able to achieve consistent precision of less 
than 1 mmag. 

By combining the single-channel PMT techniques with a high-performance CCD, we have rou- 
tinely achieved submillimagnitude relative photometry with respect to one or more comparison stars 
under photometric conditions. Unlike with PMTs, we are able to do simultaneous sky measurement 
and can take advantage of high quantum efficiencies. We slew between multiple stars of similar 
spectral types and apparent magnitude, but separated by several degrees on the sky. Our snap- 
shot method combines (1) semiautomated observations while rapidly slewing between target and 
comparison stars, (2) telescope defocusing to achieve the requisite photon counts, and (3) precision 
CCD photometry practices (discussed later). Here, we present data obtained using the SuperNova 
Integral Field Spectrograph (SNIPS) on the University of Hawaii UH2.2 m telescope at Mauna Kea 
Observatory. In §[2] we describe our observations, including the instrument, telescope, and targets. 
In § |3] we quantify the major sources of error, and in § |3]we discuss our efforts to mitigate these 
sources of noise. In § [5] we combine our noise sources using a simple noise model and compare with 
the noise from our observations. In § [6] we discuss our findings in the form of guidelines for other 
observers, and we conclude with a simulated example of the technique's intended application; a 
search for shallow (>2 mmag) planetary transits around stars that exhibit a Doppler signal. 
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Observations 



SNIFS is mounted on the south bent Cassegram port of the UH 2.2 m telescope located atop 



Mauna Kea. SNIFS is an integral field spectrograph (jAldering et al 



2002 



Lantz et al.ll2004l ) orig- 



inally designed for spe ctrophotometric ob servations of Type la supernovae for the Nearby Super- 
novae Factory project (jPereira et alJl2010l ). SNIFS consists of blue and red spectrograph channels, 
along with an imaging channel, mounted behind a common shutter. The imaging channel consists 
of two E2V 2048 x 4096 CCD detectors (CCD44-82-B23), one a high-quality science-grade device 
dedicated to imaging and the other an excellent engineering-grade device used for guiding. Each 15 
lim pixel subtends 0.137" and the field of view is 9.35'on a side. The gain for the relevant amplifier 
of the imaging detector is 1.58 e~/DN and the readout noise is 4 e~, or 6 e~ if electronic pickup 
noise is included. The telescope control system is accessed by the SNIFS instrument computer, 
allowing the observer to run scripts that i nclude telescope slews and offsets for target acquisition, 
focus changes, and guiding as described in lAntilogus et al.l (|2008l ). 



Our target and comparison stars are bright {V < 13) late K and early M main-sequen ce stars 



that are the targets of a combined Doppler-transit planet search (jLepine and Gaidosll201ll ). These 



stars are drawn from the proper-mo tion-selected SUPERBLINK catalog survey based on V-J color 
and color-magnitude relationships (jLepind l2005l ) . Candidate targets are first screened using low- 
resolution spectra to confirm spectral type, det ermine metallicity, an d remove those with strong 
Ha emission (potentially active or flaring stars, iKowalski et al.ll2009l ). Stars are the n monitored 



for b arycenter motion with the High Resolution Spectrometer on the Keck I telescope (jApps et al 



2010l ). Targets showing significant radial velocity variation are then monitored for transits. The full 
catalog comprises more than 13,000 M and K dwarf stars, and thus the nearest available comparison 
star from the catalog is usually no more than 2-3° away. 

Observations were obtained between 2010 June 22 and 2011 January 15, usually within 7 days 
of a full Moon. Photometric observations presented here usually took up approximately half of 
each observing night. Observing conditions are presented in Tabled) Seeing was measured each 
time the telescope was focused: it is the full width at half-maximum of the PSF at optimal focus, 
interpolated between discrete focus settings. Data on atmospheric extinction were obtained from 
the SkyProbe of the Canada- France-Hawaii Telescope (CFHT), which images Tycho stars around 
the boresight of that telescope through B and V filters, and d etermines extinction by comparing 
the observed signals to reference values (jSteinbring et al.ll2009l ). We report the median, standard 
deviation, and 95 percentile values of all the measurements in a night, excluding the first 10 and last 
10 observations made near sunset and sunrise. The SkyProbe data often contain high-extinction 
artifacts created when the telescope slews, so we check the SkyProbe data against videos from 
the UH 2.2m all-sky camera^, which can provide visual confirmation of clouds. While extinction 
during the nights of June 22-24, July 27-29, and January 15 was low and constant (extinction of <0.2 
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mag is considered photometric), extinction during June 27-30 and January 15 was more variable; 
variation in extinction is more detrimental to relative photometry than the absolute value of the 
extinction. Only June 30 showed evidence of (thin) clouds in the all-sky camera, predominately 
during the first few hours of the night. We also report the estimated precipitable water vapor 
(PWV) above Mauna Kea Observatories for each night, based on the mean of the 225 GHz optical 
depth {T225) measurements from the Caltech Subni illimeter Observatory, and using the conversion 



PWV= 20 X (r225 — 0.016)mm (jPavis et al.lll997l ). These values are germane because of H2O 



absorption features between 6000 and 10,000 A, i.e., in the SDSS r, i, and z passbands. . 

A single observation cycle consisted of (1) the acquisition of a smaller (800 x 800) image to 
locate the star after a slew, (2) calculation of the centroid of the focused stellar PSF, and the 
required offset to place the star at the desired location on the CCD, while at the same time (3) 
defocusing of the telescope by a specified amount, (4) acquisition of the image for photometry, (5) 
slewing to the next target while the image is read out, and (6) refocusing the telescope. To minimize 
readout time, only the bottom 1024 rows of each CCD/amplifier combination were used to image 
the target. For short slews, a single observation cycle requires ~2 minutes to complete. Nearly all 
observat ions were performed through the SDSS z passband, although some were performed with 



SDSS r (|Fukugita et al.lll996l ). 



Observations are done in one of four different semiautomated modes, (1) a target with a 
comparison star(s) in the same field of view, (2) single-set snapshot, consisting of a single set of 
stars (usually one target and one to two comparison stars) of similar spectral type and magnitude, 
and permitting the shortest cadence, (3) survey snapshot, where we observe two or more sets, with 
the stars in a set separated by less than a few degrees, but with the sets separated by as much as 
30°, and (4) constellation mode, where we repeat a sequence of observations of a large number of 
neighboring stars in a circuit, each star serves as a potential comparison star for every other star 
in the constellation (Fig. [2|). 



3. Sources of Noise 

When the target and comparison stars are in the same field of view, differential photometry 
between the two is performed through nearly the same air column. In snapshot photometry of stars 
separated by a degree or more, observations are neither simultaneous nor made through the same 
path in the atmosphere. Inhomogeneities and atmospheric transmission fluctuations (ATFs) will 
produce errors, but these can be mitigated by judicious choices of pas sband, allowable ai r mass, 



and comparison star, as well as information about the atmosphere itself (jStubbs et al.ll2007l ). Here, 
we discuss each noise source, including an estimate of the amount of noise from each source in our 
observations. We discuss some basic strategies for minimizing these sources of error in § U and in 
§[5] we combine all noise sources to predict the noise level for each observation to be compared with 
the actual precision. 
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To compare with snapshot observations, we twice conducted an experiment in which two 
stars of similar spectral type and z-band magnitude that fortuitously fell within the same field 
of view were continuously observed (cadence, 2.4, 5.5 minutes). The rms of the photometry was 
9.7 X 10~^ and 8.4 x 10~^ (Table [2]). Noise from these observations is primarily due to scintillation, 
noise from the detector (see ^ ^3.11 and 13.41 respectively), and Poisson noise. In these experiments 
the comparison stars were of comparable fluxes as the target (total counts ^■IxlO^ e~). Total 
signal detected in both tests was similar to that of our snapshot observations, where median signal 
is ~1.4xl0^ e~ for both target and comparison stars, yielding Poisson noise of ~2.7xl0~^ for 
each integration. The precision for these experiments was not better than that for our snapshot 
experiments. 



3.1. Scintillation Noise 



Scintillation noise from a telescope on Mauna Kea (4200 m elevation) is approximately 



5500 



mmag. 



(1) 



where D is the diameter of the telescope in meters, z is the air mass, t is the inte gration time 
in seconds, and A is the wavelength in angstroms (lYoungi Il993l : iBirney et al.l 120061 ). Equation 
([l]) is onl y approximate, as scintillation a. l so varies seasonally , with the speed and direction of 



the wind (lYoung 



1974: 



Dravins et al. 



1998 



Birney et al.ll2006l ) and with conditions in the upper 



atmosphere (jHeasley et al.lll996l ). For example, the true value of the air mass exponent is ~2 when 
observing in the same direction as the wind and is ~1.5 when observin g perpendicular to it, and the 
wavelength dependency vanishes for la rger telescopes (IRoddieii 119811). M o re rig orous methods of 



calculating scintillation can be found in lKenvon et al.l (j2006l ) and iKornilovl (j201lh . but with proper 
modification of the air mass term for wind direction and no wavelength de pendency, equation 



(fTl) is consistent within 20% of scintillation measurements from Mauna Kea ( Daintv et al. 



Gilliland et al 



1993 



1982 



Cherubini et al.ll2008l ) under conditions similar to those of our observations. 
With wind modifications and without the wavelength term, equation ([T]) gives a median value 
per star of 2.4 X 10""^ for noise due to scintillation in our observations. It becomes the dominant 
contribution to the error budget at exposure times if less than 15 sec. 



3.2. Extinction 

First-order extinction: Our target and comparison stars, although close on the sky, will have 
a nonnegligible difference in air mass. For an observation of stars lasting several hours, the differ- 
ential air mass of the target and comparison star(s), x = Xtarg — ^comp, will change. For small 
values of X; and assuming no significant change in the extinction coefficient, the fractional change 
in the normalized relative fiux will be 6f ~ 0.4 ln(10)ii^(A)Ax, where E{X) is the extinction co- 
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efficient in magnitudes and air mass and Ax is the change in air mass difference over the course 
of an observation. A typical value of A^; in 2 hr for a well-selected comparison star is 0.003. As- 
suming an extinction coefficient for photomet ric conditions of ~0.04 mag air mass~^ in SDSS z 
(jMcCord and Clarklll979l : iHodgkin et al.ll2009l ). the change in flux of the target star relative to the 
comparison star will be ~10~^. Even for comparison stars within 5° of the target, values of greater 
than 0.01 for Ax are possible, and the resulting change in flux could be as much as ~4xl0~'^. If 
large values of Ax combine with higher extinction coefficients (i.e., nonphotometric or high mean 
PWV), the systematic effect can be larger than ~10~^. However, the resulting trend is systematic 
and can be removed with ~10~'^ precision provided that the extinction coefficient is measured (or 
known) to better than 0.01 mag air mass~^. 

Second-order extinction: If the target and comparison stars do not have the same spectral 
energy distributions over the passband AA, changes in the extinction by the atmosphere with 
air mass and time will produce trends in their relative signals. The ef fect will be proport ional 
to AA^. Second-order ext in ction error was ana lyzed by (among others) I Young et al.l (jl99ll ) and 
Everett and Howelll (120011 ). lYoung et al.l (Il99ll ) estimate the size of the effect by observing two 
stars differing in B-V by 0.3 mag and find that the difference in B extincti on of the two stars 
i s ~7 mmag air mass~^, although other estimates are notably smaller: e.g., lEverett and Howell 



(j200ll ). The color effect is smaller in redder passbands, with stars of later spectral type, and with 
narrower bandwidth filters. In the case in which the spectral energy distributions of the target and 
comparison stars are known, a correction can be made. 

We estimate second-order extinction errors with respect t o a typical tar get star in our observa- 
tions (4000 K and log g = 4.5). We convolve rnodel s pectra of lKurucj (jl99ll ) with the transmission 
curve of the SDSS r, SDSS z (jFukugita et al.lll996l ). B, or a narrow custom z (J. Johnson 2011 
private communication) filter, the quantum efficiency of the EV2 detector, and a model of atmo- 
spheric transmission. Fig. [3] shows the profile of each of these filters multiplied by the SNIFS 
transmission and with an approximation of the atmosphere above Mauna Kea. Our atmosphere 
model is based on l ow-resolution spectra of standard stars taken with SNIFS, combined with the 
HITRAN software (jRothman et al.ll2009l ) in regions where our spectra are incomplete or contam- 
inated by features of our standards. Fig. [4] shows the expected change (systematic error) in flux 
ratio {Ftarget/ Fcomparison-, normalized to 1), as a function of air mass, for the SDSS r, SDSS z, 
B, and narrow z filters. We assume that the target and comparison stars are observed through 
identical air mass. We find the size of the effect in the B band to be approximately 2.5 mmag 
air mass~^ for a c omparison star with (B-V) color difference of 0.3, which is significantly smaller 
than that found bv I Young et al.l (|l99ll ). The difference is likely due to the difference in atmosphere 
above the observatories (Mauna Kea versus Mount Hopkins), different instrument proffies, and the 
choice of later spectral type stars for our calculations. For our observations in the SDSS z filter, 
we calculate a median value of only 4.5 x 10~^, which is a benefit of choosing comparison stars of 
a similar spectral type to the target star and working at longer wavelengths, where the differences 
in stellar spectra are smaller. 
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3.3. Short-Term Atmospheric Transparency Variations 



Snapshot obse rvations are not simultaneous (> 2 min lag), and variations in the atmosphere 
will produce error. iHill et al.l () 19941 ) find the average power spectrum of transparency fluctuations 
in the atmosphere above Mauna Kea to obey 



logP(z^) = -9.84- 1.50 log z/, 



(2) 



where v is the frequency in hertz. Equation Q is based on an average over 691 clear days. An 
individual night could have significantly higher or lower ATFs, which makes this source of noise 
particularly difficult to estimate for any given observation. Regardless, the spectrum contains no 
coherence time within the timescale of our observations, and the photometric stability improves 
with decreasing time between target and comparison star observations (shorter cadence), although 
with diminishing returns. We make use of equation ([2]) to estimate the noise from atmospheric 
variations between observations of target and comparison stars. In general, the variance due to 
ATF is 



(7 ATP 



D{t)S{t)dt 



(3) 



where D{t) and S{t) are functions that describe the transparency fluctuations and integration 
windows respectively, and <> represents the time-averaged expectation. We assume that our 
observations behave as a function of the form 



Sit) = I 



1 







t £ Tcomp 



(4) 



where Ttarg and Tcomp are the times over which we are integrating on the target and comparison star 
respectively. We can represent S{t) and D{t) as Fourier transforms ^(w) and ^{oo) and substitute 
into equation [3l After some simpliflcation, we flnd 




{&{u)f{y{u)fdu 



(5) 



where i^(a;)^ = P{io), and uj' is the angular frequency over which &{uj) and ^(io) are defined. 
Equation ([2]) covers a significant portion of the optical spectrum, whereas we are interested in the 
z band, where th e atmosph e re is more transparent. We assume that P{u}) at the z-band scales 
linearly with the iHill et al.l (jl994l ) power spectrum and we derive a scalar coefficient based on 
our atmosphere model discussed previously and information gathered from the CFHT SkyProbe 
(jSteinbring et al.ll2009l ). Our atmosphere model gi ves us the averag e atmospheric transparency for 
SDSS z (or any other bandpass) to compare with iHill et al.l (jl994l ). We use the CFH T SkyProbe 
data to derive a scale factor between transparency fiuctuations of iHill et al.l (|l994l ) and that of 
the y-band. This is consistent with the results of our atmospheric model; i.e., the amplitude 
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of transparency fluctuations scales with average transparency fluctuations, thus enabling us to 
properly adjust equation ([2]) to any bandpass within the range of our atmosphere model. For our 
observations, the median noise due to transparency fluctuations is 2.8 x 10"'^ and is typically smaller 
than Poisson and scintillation noise. 



3.4. Detector Noise 

No CCD has a perfectly uniform response rate, even when corrected with flat-fields. We 
calculated the error due to pixel response nonuniformities combined with motion of the defocused 
image around a given coordinate centroid, ignoring changes within the defocused PSF and chip 
nonlinearities. To accomplish this, we took 225 dome flats in 2010 June with exposure times from 
2 to 15 s; 150 flats were taken using the SDSS z filter, and the remaining 75 were taken with 
SDSS r. All flat-fields were obtained within a period of 5 hr. Count levels varied from 5 x 10^ 
to 3.8 X 10^ e~ (typical levels for the signal in one of our defocused images). For each filter, we 
performed a linear least-squares fit of electron counts in each pixel versus median electron counts 
(for the chip of interest), assuming Poisson and read-noise variance. The result was a pixel- by-pixel 
map of the CCD's response for each filter. The median count was used instead of exposure time 
to remove effects of lamp variation. We created a model defocused PSF and scanned it over the 
1024^ region of interest in our detector response map to find the region(s) of the CCD with the 
most uniform pixel response. In general, the CCD's behavior was similar in both filters, although 
variation of the pixel response for the best (most uniform) regions was slightly better when using 
the r filter (rms pixel response ~ 0.7%) than it was when using the z filter (rms ~ 0.8%). For 
either filter, even small PSF motions (< 3 pixels) on a bad region (very nonuniform pixel response) 
can cause noticeable (> 3 x 10~^) variations in received flux. Large motions (> 15 pixels) around 
a well-behaved region (highly uniform pixel response) can create similarly large variations in flux. 
Because we use a good region on the chip, typical rms motions for our observations (< 5 pixels) 
contribute minimal noise (^ 1 x 10^"^ for z and < 0.8 x 10~^ for r) to the total error budget (Fig. [5]). 

Pixel-to-pixel variations in the response of the detector can be removed to a certain extent 
with dome and/or sky flats. However, using flats introduces additional noise to the data, because 
quantum efficiency variations among the CCD pixels are wavelength-dependent, the dome lamp 
(or sky) will not have the same spectrum as the target, and master flats composed of numerous, 
high-S/N flats still have noticeable Poisson noise. Twilight flats will better match the spectral 
energy background distribution of the data frames; however, it is difficult to get a large number 
of high-S/N twilight flats in the relatively short twilight window. These and many othe r error s 



associ a ted w i th fiat-fiel d ing a re dis cussed more thorough ly in (among others) iNewberryl ()l99ll ) 



TobinI (ll993l ). lManfroid (|l995l ^. and lManfroid et al.l (|200ll l. 



Noises brought about by fiat-fielding may be significant, compared with the < 10~^ sized errors 
induced from inaccurate detector response (Fig. [5]). For example, consider a master flat composed 
of 10 twilight flats, each with ~ 20000 e^ pixel~^, and defocused target/comparison stars spread 



- 10 - 



over 1000 pixels. Taking into account Poisson noise only, the error associated with flat-fielding to 
the light curve will be ~ 1 x 10""^, putting it on par with other noise from detector nonuniformity. 
Further, when we apply a median of 10 twilight flats to our dense grid of 150 flats (we assume the 
grid of flats to be a more accurate map of the detector) the predicted errors created by image motion 
do not improve signiflcantly, as shown in Fig. [5j Applying a flat-field correction improves precision 
when using a mediocre region of the chip, and/or when the image is drifting ^10 pixels. Since we 
make use of a very flat region of the chip, and image motions are small, flat-field corrections may 
actually add noise. 

Detector nonlinearities are small for most science-grade instruments, but they can become 
important for photometry at submillimagnitude precision. For a single observation the recorded 
fiux will be approximately 



where Xi is the normalized incident intensity on any given pixel, and a is the level of nonlinearity 
in the pertinent range. If the incident fiux is x + r, where r is normally distributed with standard 
deviation a, then the difference between the signal of any two observations, Sj and Sk, will be 



The distribution of Sj — will have an average of zero, but a nonzero standard deviation. If we 
ignore terms of higher order than cr^, then the latter will be: 



where N is the total number of pixels. A typical aperture size for our observations contains 
N ~1.5xl0^ pixels, a can vary signiflcantly between observations, but is no higher than 1.0 
and is often less than 0.1. In 2005 a series of dome flats were taken to determine the flux range 
over which the SNIFS E2V CCD is linear, which we use to estimate a. Exposure times varied 
between 2 and 55 s and were interleaved to reduce any time-varying effects. All exposures were 
taken within a few hours. The region containing the lower 200 pixels of the photometric chip and 
relevant amplifier was found to be linear to 0.15% ± 0.08% between 7 x 10^ and 1x10^ e". The 
rms of these measurements is 0.16%, most of which is due to variations in lamp brightness. Since 
we require the chip to be linear to better than 0.1% for submillimagnitude photometry, we aim to 
keep the maximum fiux in any given pixel below 5.0 x 10^ e~, where the CCD is linear to better 
than 10~^. Thus for fiux levels consistent with the majority of our observations, a < 0.0015. Using 
1.0 as an upper limit on a, the noise due to nonlinearities is only 1.7 x 10^^. 

Another potential source of error is variation in the shape of the PSF, combined with fixed 
calibration errors in the pixel response. Even with sufficient defocusing, phase errors from seeing 
will induce fluctuations across the face of the image. The result is that, although total fiux may 
be conserved, fiux will be redistributed between pixels, which will produce uncertainties to the 




(6) 




(7) 



rms ~ 



^/2acj2 



(8) 



- 11 - 



extent that pixel responses are not perfectly calibrated. If the PSF did not vary, such calibration 
uncertainties would not produce time variations. Fig. [6] shows the rms difference between eight 
individual images and the mean. The stellar wings are mostly constant, but there is significant 
variability in the shape of the core annulus, particularly at its inner and outer edges. We estimate 
the magnitude of such an error by assuming that the calibration suffers from errors having a fixed 
Gaussian distribution with a = 1%. The error calculated from the eight images is only 4 x 10^^, 
which is a consequence of summing uncorrelated errors over the more than 10^ pixels with which 
the signal is acquired. 



4. Strategies to Minimize Noise 

4.1. CCD Regions and Signal Levels 

Large pixel nonlinearities, like those seen in some near-infrared detectors, can be corrected 
by applying fits t o an intensity seri es for each pixel and then implementing the corrections to 



the detected flux (jVacca et al.l |2004| ) . Pixel nonlinearities can also be mitigated in differential 
photometry by using the same pixels and similar flux levels and exposure times. More subtle 
nonlinearities and time- variable pixel responses are difficult to independently measure and remove, 
but cannot be ignored when submillimagnitude precision is required. One solution is to identify 
those pixels that are poorly behaved (highly nonlinear and/or time- variable pixels) and avoid them. 
It is not possible to simply mask these pixels during image processing, because that sensitizes the 
total signal to small changes in the position or shape of the PSF. Our strategy for dealing with 
noise related to the detector is threefold, we locate the range of signal over which the detector is 
sufficiently linear (discussed in § 13. 4p , identify regions of the detector containing pixels with uniform 
response rates, and minimize the motion of the defocused image between integrations. 

We used the series of flats taken in 2010 (discussed in § 13. 4p to locate time- variable or highly 
nonlinear pixels on the SNIFS E2V detector. We identified pixels that have the highest values 
for a linear fit. Fig. [7] shows a histogram of the values for pixel across the 1024^ region of the 
detector used for photometry. Because our flat experiment is designed to locate the best and worst 
regions on the detector the absolute variation from linearity is not as important as the identification 
of pixels that show higher levels of nonlinearity (higher x^) than others. ^0.1% of pixels have 
values above the expected distribution (x^ > 240) and are predominantely concentrated in a few 
areas (e.g., near the edges of the detector), making poorly behaved pixels easier to avoid. 

Our data indicate that the best precision is achieved without flat-fielding, provided that images 
are placed in good (flat) regions of the chip and that image motion is small (rms < 10 pixels). 
Fig. [8] shows a comparison of precision with and without flat-fielding. For light curves built with 
flat-fielding we use twilight flats where available, and dome flats elsewhere. Although some light 
curves benefit from the use of flat-fields, applying fiat-field corrections to most of our data worsened 
the precision. Instead, noise from errors in pixel response are minimized by using the same (good) 
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region of the chip for all observations with a given filter (see Fig. [5] for a comparison of noise 
from good and bad regions of the CCD) and by keeping image motion small from integration to 
integration. SNIFS' automated acquisition and position of the target/comparison stars is more 
than sufficient to keep image motion below 5 pixels. Flat-fielding will have more beneficial results 
for detectors with higher interpixel QE variations or higher image drift between integrations. 



4.2. Choice of Filter 

Choice of bandpass filter is a tradeoff between the desire to maximize signal from the star 
and to minimize the contribution from wavelengths that are affected by scintilla tion, scattering 



(by aerosols and molecules), and absorption by O2 and H2O (jStubbs et al.l 120071 ). Observations 
through redder filters will have lower first- and second-order extinction (see ^3.2p . Stellar variability 
i s a significant source of noise in the V b and for as much as 20% of M dwarfs on 12 hr timescales 



(jHartman et al.l 120111 : iCiardi et al.l 120111 ). However, observing at longer wavelengths reduces the 
spot-photosphere contrast, reducing noise from stellar rotation and spots. Among the redder SDSS 
filters, i is seriously compromised by molecular bands (Fig. [3]), 2: is bracketed by two H2O bands, 
although the low quantum efficiency of most CCDs redward of 9000 A mitigates H2O contamination, 
while r is more affected by aerosol fluctuations. Use of a narrow z filter (if one is available) would 
almost completely eliminate second-order extinction error (Fig. H]), due to its narrow bandwidth 
and red central wavelength (Fig. [3]). A narrow z filter would also reduce noise from scintillation 
and scattering from aerosols while avoiding a major H2O line. 

Observations in near-infrared filters are more affected by atmospheric emission than bluer 
bands, however, our data suggest that the overall noise from atmospheric emission and/or fringing 
is small. Fringing is present in SDSS i and z for SNIFS, but noise from fringing is minor compared 
with the rest of the sky background. Since we predominantly observe within 7 days of a full Moon, 
the dominant component of the sky background is scattered light from the Moon. Fig. [3] includes 
an estimate of the atmospheric emission over Mauna Kea based on spectra taken from SNIFS. In 
the z band, integrated sky emission in a ~ 7 " radius, a typical defocus size, will be at least a factor 
of 100 smaller than the fiux from a typical star (m^ = 9) in our program. In the narrow z band, 
this grows to 180, since the narrow z-filter transmission cuts off before the largest OH emission 
lines that contaminate the z band. 

To test the performance of SDSS z and r filters we performed regular observations while 
interleaving exposures of r and z filters. We kept other controllable variables (exposure time, 
defocus, etc.) fixed between filter changes. The z band slightly outperformed the r filter overall 
(see Table [2|), but not by as much as our total noise estimate predicts (see §[5] for more on the noise 
model). The discrepancy could be due to underestimating the effect of the H2O line on the z filter, 
which is dependent on the PWV of that particular night and induces noise that cannot be easily 
modeled. 
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4.3. Degree of Defocusing 

A minimum degree of defocusing is required for a fixed total signal from the star and linear 
range of the detector. For example, if the linear range of a CCD is 5 x 10^ e~, then the acquisition 
of 2 X 10^ e~ must take place over at least 400 pixels (for SNIFS, this is a circle 5" across). Because 
the PSF is nonuniform, more pixels must actually be used to collect the signal and avoid saturation. 
Greater defocusing may be desirable for very bright stars to avoid integration times much shorter 
than the CCD read time and to minimize scintillation noise. 

Fig. [6] is the mean of eight 30 s integrations of a defocused y ~ 13 M star obtained through 
the SDSS z filter. Most of the signal is confined to an annulus with a ~4" radius that represents 
an out-of-focus image of the telescope's primary mirror. The defocused PSF core is surrounded 
by a halo pattern, produced by the defocused convolution of the telescope pupil with the seeing, 
which extends out ~10". The image is i iot axisymm etric because of coma in the telescope optics 



(C. Aspin, private communication 2010: lBehiill973l ). The area within the annulus is ~10 pixels 



and the average signal per pixel is 1.7 x lO'^ e , with a total signal of 1.6x10^ e . The total read 
and electronic pickup noise (200 e~) is much less than the photon noise (4x lO^e"). The maximum 
S/N is achieved by summing the signal within an aperture of 8"radius, if only photon noise, sky 
noise, read noise, and electronic pickup noise are considered (Fig. [9]). 

More defocused images are less sensitive to motion of the stellar image on the CCD; however, 
the total read noise and probability of including ill-behaved pixels are greater. We consider the 
idealized case where the PSF is a circular aperture comprising pixels, each of which gets S x N^^ 
signal (e~) and has read and electronic pickup noise a\ (e~ pixel"^) and uncorrected response 
variation (T2, which we presume adds in quadrature. The variance from a single nonover lapping 
pixel due to the last source of error is S'^N~'^a\. If the centroid rms motion in pixels is b then the 
number of pixels that are not common to a pair of stellar images is '^2it5^J N /tt. The total variance 
due to read noise and image jitter is 

^2 ~ Nal + 2A-3/2^55Vi. (9) 

Total read noise increases with A, but noise due to image jitter decreases. The number of pixels 
A^, that minimizes the total error is 

N^^{25alS^a^^f'\ (10) 

For 5 = 3 (the median centroid motion in pixels for our observations), a2 = 0.01, ai = 6e~, and 
S = 10 ''e~, and A,^ ~ 7.2 x 10^: i.e., a circular region with diameter of ~96 pixels (13 ") with 
total noise a = 6.6 x 10^^ e^. A* is sensitive to well depth, read noise, and the uncorrected pixel 
response noise, but can be estimated by experiments with defocused images. We are often forced 
to use higher defocus values than equation ([TU]) suggests to keep the counts in the coma-induced 
peak in the linear range of the detector. 
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4.4. Integration time 



Poisson and scintillation noise decrease with integration time whereas noise from ATF will 
increase. Additionally, higher exposure times for a given star will require more defocusing, increas- 
ing contributions from read and sky noise. There exists an optimal exposure time that minimizes 
the total error from these sources. We estimate the optimal exposure time as a function of the 
magnitude of the star using our calculations for all of these sources of noise: ATF, scintillation, 
Poisson, read and electronic pickup, and sky. We use a model of our defocused images (including 
the coma-induced peak) to calculate the required defocus in order to satisfy equation ([TO]) and to 
keep the flux below 5 x 10^ e~/pixel (i.e., where the detector is linear to better than 0.001). We 
then calculate the total expected noise as a function of exposure time, assuming target and refer- 
ence stars of identical flux and spectral type, yielding an optimal exposure time for a given stellar 
magnitude. We repeat this calculation for a range of magnitudes encompassing those in our survey 
(see § [2|). Fig. [TO] shows the optimal exposure times, defocused image radius, and calculated noise 
for the exposure time with the best expected precision. Our calculation overestimates exposure 
times for the brightest stars (m^ < 8), where the suggested level of defocusing becomes impractical. 



In our exposure time tests, we interleave observations of a pair of stars {niz = 10.8 and 11.4) 
using three sets of exposure times each (25s, 50s, 100s and 40s, 80s, 160 s, respectively). According 
to our applied noise model, the best precision should be obtained with the 50s and 80s integration 
times (Fig. [T0|) . However, the lowest total noise came from the highest exposure time (100s and 
160s). Additionally, the 50/80s and the 100/160s exposure time sets beat the expected precision 
(see Table [2]), suggesting that ATF may have been lower on that particular night, or that we 



might be averaging oyer stellar P-mode oscillations with lo nger exposure times (| Young et al 
Brown and GillilandlflQoi : lHenrvlll99d : Iciardi et allboill ). 



1991 



4.5. Choice of Comparison Star 

There are noticeable gains from proper selection of a comparison star. Our collection of low- 
resolution spectra enables us to select comparison stars of a similar spectral type to our target, 
resulting in very small errors from second-order extinction (4.5 x 10^^). If we replace all comparison 
stars with G dwarfs (6500 K), the median value becomes 3.0 x 10~^, which would make it a more 
significant part of the error budget (Table [3]). This is especially important with M-type target 
stars, which are intrinsically faint and are therefore less likely to be selected as comparison stars 
based on flux alone. 

In principle, more comparison star observations will improve the S/N, but these additional ob- 
servations will occur at increasingly earlier or later times and can engender greater systematic error 
due to changes in the atmosphere. We examine the constellation approach by linear construction 
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of a reference signal s for the jth observation of the ith star: 

Sij = "^^dijklSkh (11) 
kl 

where Ski is the flux, and the weighting function aiju = for i = A:. In general, for nonzero values, 
aijki will depend on the characteristics of the atmosphere, the cadence of the observations, and the 
stability of the data from the stars (i ^ k) in the constellation. To keep the signal normalized, for 
all i, j we force: 

'^aijki = l- (12) 

kl 

We use the constellation technique that produces the best 

photometric precision. We compute a grid of reference signals for each star following equation (|lip 
and apply a range of possible values for ajj^/. For our observations, we find that the best precision is 
achieved most consistently when using a reference signal formed from the stars immediately before 
and after the target observation, i.e., a^j^i ~ for all \i—k\ ^ 1 or j ^ I. This conclusion is consistent 
with our estimates of the noise from ATF (§ [33]), which begin to overwhelm the noise budget when 
more than ~ 7 minutes passed between target and comparison star observation. However, the 
additional stars in the constellation can be used to identify and mitigate noise from aberrant data 
points (e.g., data contaminated by cosmic rays), or variable stars in the constellation. Fig. II II shows 
the precision for light curves built using the stars before and after the target observation alongside 
light curves (of the same stars) constructed using an equal weight of all available comparison stars. 



5. Comparison to Observations 
5.1. Total Errors 

We calculate the total predicted error of each observation by adding the uncorrelated noise from 
both the target and comparison star(s) due to: scintillation, Poisson (of the star), sky (background), 
motion on the chip, change in shape of the defocused image, overall chip/amplifier nonlinearities, 
and readout (read plus electronic pickup noise), along with first- and second-order extinction and 
ATF. Noise from motion on the chip is estimated using our model of detector behavior (§ 13. 4p . 
The noise from the change in shape (changes in distribution of flux over the defocused image) is 
calculated assuming that the detector interpixel response suffers from 1% Gaussian variation. We 
correct for noise from first-order extinction when this noise source is greater than 1 x 10~'* and we 
know the extinction coefficient to better than 0.01 mag air mass^^. Although specific levels for each 
term vary significantly between observations, the median values of Poisson noise are the largest, 
followed by noise from ATF, scintillation, first-order extinction, sky, read and electronic pickup, 
motion on the chip, second-order extinction, changes in shape of the PSF, and overall/amplifier 
nonlinearities (typical noise levels are listed in Table [3|) . The small relative size of many of these 
errors (second-order extinction, for example) is mostly due to judicious choices in each of our 
observations: e.g., wise choice of comparison stars, exposure times, and region of the detector used. 
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Fig. [12] shows the precision of each observation versus the total theoretical noise for each 
observation. Although the mean of the predicted error is within 3% of the theoretical precision, 
there is significant scatter (rms = 12%). We expect the major discrepancy to be in our estimate 
of error from ATF, which is difficult to calculate for any specific observation. Other noise terms 
are approximate, such as scintillation, which is only good to ~20% (even with modification for 
wind direction). We have also assumed no correlation between terms, which might be causing 
us to systematically underestimate noise. Additionally, stellar variability may contribute to the 
noise budget, e.g., P-mode oscillations, whic h can be noticeable for our targets (mostly late K 



and early M dwarfs) at our level of precision ([Brown and Gillilandlll994l : [Henrv[[l999l : [Ciardi et al, 



2011 



). Observations with short exposure times are especially susceptible to P-m ode oscillations, as 

Thus, our model 



P-mode oscillations have timescales on the order of minutes ([Young et al. 



1991 



correctly accounts for the overall magnitude of errors but cannot reliably predict the error of a 
specific observation. 



5.2. Correlated (Red) Noise 



The performance of most transit surveys is limit ed largely by tinie-correlated (red) noise, which 
limits gains in precision from binning/phasing data ([Pont et al.[[2006l : [Carter and Winn[[2009l ) . Red 
noise is also especially deleterious because it can mimic a signal of interest, leading to precise but 
inaccurate results. We expect to have some degree of time correlation in noise from extinction, 
ATF, and scintillation, which will manifest as a combination of white and red noise. Due to all 
these considerations, we analyze the level of red noise in our data. We estimate red noise from 
ou r comparison-star-correc ted light curves using a wavelet-based method. The method, described 



m 



Carter and Wind ([20091 ). determines parameters of noise formed as a combination of Gaussian 
(white) noise and noise with a power spectral density varying as l/f^. We test the simple case of 
7 = 1, as well as the case of 7 = 1.5, implied by equation ([2]). No light curve shows red noise above 
6 X 10~^ for either 7 = 1 or 7 = 1.5, and median red noise is only 2.8 x 10~^, at least some of which 
is due to stellar variability. The low levels of red noise compared with total noise suggest that 
we are not being overly hindered by time-correlated noise sources. Thus data acquired using the 
snapshot technique can be binned to produce improved precision; after binning each observation 
to a cadence of less than 20 minutes, the median rms from our observations is 3.2 x 10~^ (all over 
intervals of at least 2 hr to ensure sufficient data points) and all binned light curves had rms of 
less than 10"'^. Further precision improvements can be made with additional binning; however, a 
cadence of greater than 20 minutes is impractical for detecting the presence of a transit. 



6. Discussion 



Using defocused snapshot photometry, we have consistently achieved submillimagnitude pho- 
tometry from the summit of Mauna Kea under clear conditions. Our lowest photometric rms was 
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5.2 X 10~^ (with a 5 minute cadence) over a ~2 hr interval. Of our 38 experiments using the 
snapshot method, 32 of them had rms precision of less than lxlO~^, and our median precision was 
7.8 X 10~^. We include our constellatio- mode observations, where we observe eight stars in se- 
quence (Fig. [TTI). Al t hough our precision i s not as good as the best recorded from the ground (e.g., 



Gilliland et al.l Il993l : iHartman et al.l |2005| ) , it is sufficient to detect the transit of a Neptune size 



planet around a dwarf MO star or a lOAf^ super-Earth around an M4 star at 5a significance. Fur- 
ther, our technique does not rely on a large number of comparison stars, yields submillimagnitude 
photometric precision consistently, and can be readily automated. 



6.1. Methods of Submillimagnitude Photometry 

Although our methods were optimized for the UH 2.2m telescope, SNIFS instrument, and 
Mauna Kea observing site, we give guidelines for achieving submillimagnitude pho tometry that 



can b e foll owed by any one. Some of these are recapitulated guidelines offered by lYoung et al 



(| 199 il l and iHenrvl (|l999l ). 



Appropriate selection of comparison star: Optimal choice of comparison star is a tradeoff be- 
tween minimizing the separation on the sky and using a star of comparable apparent brightness and 
spectral type. Submillimagnitude photometry becomes significantly more difficult with separations 
of more than a few minutes of slew time between target and comparison stars. Using more than 
two comparison stars provides no significant improvement (and usually yields inferior precision) , as 
changes in the atmosphere grow larger than other sources of noise related to the comparison star(s) 
after just 7 minutes. Using two comparison stars is preferred over a single star in order to identify 
any signal that is the result of stellar variability instead of a transit. Large separations increase 
air mass differences, and give the atmosphere time (and distance) over which to change. Although 
stars of similar apparent brightness are preferred, this technique allows the observer to use different 
exposure times and focus settings for target and comparison stars. Thus, significantly dimmer or 
brighter comparison stars can be used with only minimal loss in precision. The spectral type of a 
comparison star should be similar to that of the target, especially if the stars are observed over a 
wide range of air masses. It is not usually worthwhile to use comparison stars that require more 
than ~ 6 minutes of time between target and comparison star observation. Past this, the precision 
loss from atmospheric fluctuations overwhelms noise from a poor choice of comparison star. 

Exposure times and defocusing: Integration times and the amount of defocus must be chosen so 



that t he signal per pixel is within the linear range (<C 1% nonlinear) of the detector. lYoung et al 



(jl99ll ) suggest using shorter integrations to move rapidly between target and comparison stars. Our 
calculations and data indicate that the optimal exposure time depends on the tradeoffs between 
ATF and shot/scintillation noise, and it often exceeds 1 minute. Fig.llOlshows how to scale exposure 
times and defocus radius with the brightness of the star, although the values may be significantly 
different for other telescopes. Assuming there are no stars that will overlap with the target when 
defocused, and that the detector has a sufficiently large region of well-behaved pixels, it is best to 
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find the ideal exposure time for eacli star, and defocus as appropriate to keep flux levels below the 
nonlinear range of the detector. 

Minimize time between observations: Making use of a smaller region of the chip will reduce 
overhead from CCD readout enough that the time between observations is usually limited by 
telescope/dome slew time. In the constellation method, a number of stars are observed in sequence 
to search for transits (Fig.s [2] and llip . In order to maximize observing efficiency and the number 
of targets that can be observed with the highest cadence (> 1 observation every 30 minutes), a 
time-efficient trajectory must connect the target stars. With an efficient path, the constellation 
technique can be expanded to observe a large number of stars in sequence at the cost of cadence. If 
the dome motion is slower than the motion of the telescope, the fastest path minimizes azimuthal 
moves. 

Observe near meridian crossing: In addition to being at minimum air mass, objects near merid- 
ian crossing experience the smallest changes in air mass relative to nearby stars. When to observe 
is especially important when observing over longer periods of time (> 1 hr). For observations on 
the same set of stars lasting significantly longer than a few hours (producing large changes in the 
air mass difference), it will be necessary to take enough measurements of the extinction coefficient 
to constrain it to better than 0.01 mag air mass~^ during the observation. 

Selection of passband filter: We achieved similar photometric precision in both SDSS r and z 
passbands, possibly a result of higher scintillation noise in the former and contamination by the 
H2O lines in the latter. A narrower z filter lacking the H20-affected wings (Fig. [3]) would avoid 
both of these problems, as well as reduce the size of second-order extinction. 

Use of well-behaved regions of the detector: Regions of the detector containing nonlinearities, 
variable behavior, or significant nonuniformity must be avoided and cannot simply be corrected 
by flat-fielding (Fig. [7j). One significant advantage of the snapshot method is that images of the 
comparison star and target star can be placed in the same position on the detector. The rms of 
centroid motions (from integration to integration) should be kept to less than 10 pixels (Fig. [5]). 
Snapshot photometry only requires a field of view large enough to fit a single defocused star and 
enough surrounding field to estimate the sky background with an absolute error much less than 
that of the star. 



6.2. Performance of a Hypothetical Transit Search 

To close, we examine the performance of the photometric method described in this article 
using Monte Carlo simulations of a search for transiting planets around M dwarf stars. We assume 
a telescope like the UH 2.2m at a photometric site such as Mauna Kea. The all-sky input catalog 



their proper motions (/i > 0.15" yr ^) ( Lepine and Sharall2005l : iLepine and Gaidoal201ll ). and V-J 
colors, where J-band magnitudes come from the Two Micron All Sky Survey, a near-infrared all-sky 
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survey (jSkrutskie et al.l 120061 '). We estimate the effective temperature bas ed on an average of V-J , 
V- H, V-K, and J-K colors a nd temperature-color relations established b ylHoudashelt et al.l (12000) 
and lCasagrande et al.l (|2008l ). A bolometric correction is computed as per lHoudashelt et alj (|200Q ). 
assuming solar metallicity and surface grav ity log q = 4 .5. Absolute magnitudes are calculated 
assuming a position on the main sequence (jLepind l2005l ). We make three sepa r ate estimates o f 
the mass based on the near-infrared absolute magnitudes (jPelfosse et alJ l2000l : IXia et alJ 120081 ) 
and average them. Although a radius can be calculated from the Stefan-Boltzmann relation, we 
instead empirically infer it fr om the mass using t he relation i?* ~ M^"^ (O.SM© < Af* < O.GMq^ 



based on theoretical models (Baraffe et al 



measurements of single stars (jPemorv et al. 



IQod) and supported by the available data on radius 



20091 1. 



We use the pla net mass distribution dn/dloR Mp = 0.39M^'^^ (M i n Earths) estimated by 
Howard et al.l (120101 ) and the period distribution of ICumming et al.l (j2008l ): dn/cZ log(P)~P'''^^ for 
P = 2 — 20 00 days. These distributions are only known for solar-ma ss stars and may be different for 



M dwarfs (^Raymond et al. 



2007 



Montgomery and Laughlinll2009l ). The recently released Kepler 



catalog of planet candidates contains ~1200 objects with estimated rad ii, but <3% of them ar e 
around M dwarfs. We adopt the orbital eccentricity distribution derived bv lShen and Turned (|2008l ). 
recognizing that the distribution is w ell-constrained only for giant planets a nd that of low-mass 
planets may be significantly different (jPilat-Lohingeii l2009l : iMann et al.ll2010l ). Inclination angles 
and longitudes of periapsis are drawn from isotopic distributions. For a mass-radius relationship 
we use the surprisingly simple formula Rp = Mp'^ (Earth units), which captures the overall trend 
for exoplanets of masses intermediate to those of Earth and Saturn. This diverges from the mass- 
radius relationship of rocky planets derived from interior models (e.g.. ISeager et al.ll2007l ). possibly 
because of the tendency of more massive planets to have H/He envelopes. We assume that all 
planets above a mass of 125 Earths are the size of Jupiter (11.2 Earth radii). 

The simulated observing program is constructed as follows. Days are first selected randomly 
from throughout the year. Sunset and sunrise are determined for the Mauna Kea summit, and 
observing begins and ends 30 minutes after and before sunset and sunrise, respectively. The short 
twilight interval is justified because our targets are bright and Rayleigh scattering from the sky is 
low in the z passband. A target is selected at random from an updated subset of the catalog that 
is found at an air mass below 1.4. We allow 2 minutes for the telescope to slew to each target and 
the instrument to set up, and we assume a fixed integration time of 2 minutes. Each target that is 
observed is placed on hold and will not be reobserved until at least 2 hours has elapsed (the typical 
duration of a transit). We do not consider the requirement that each target be observed a certain 
number of times each night because we ultimately sum the probability of detecting a transit over 
all observations, not all individual targets. 

Detection requires a S/N ratio greater than 5 i n a s ingle observation. The transit signal is 
calculated using the limb-darkening model of IClareti (120001 ). The total noise consists of shot noise, 
scintillation noise, and other (atmospheric) noise (0.4 mmag) added in quadrature. The first two 
contributions are based on the performance of the SNIFS instrument and assume a SDSS z passband 
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(see § [3] for details) . The detection probability is the product of the transit probability times the 
probability of detecting the transit averaged over 30,000 Monte Carlo realizations of planets. We 
can impose a Doppler detection criterion by performing the average only over those systems where 
the radial velocity amplitude exceeds a specified value. This simulates a scenario in which targets 
are first selected by Doppler observations for transit searches. 

Fig. [13] plots the average required number of clear nights per transit discovery as a function 
of the Doppler threshold. A higher threshold selects more massive planets on closer orbits; these 
are more likely to transit and are more readily detected if they do. For a 6 m s~^ threshold, the 
mass distributio n of detected pla.nets p eaks at around 24 Me (Neptune size). If 56% of nights 



are photometric ( Steinbring et al.ll2009l ). then a 1 yr survey of 6 m systems would produce 25 
transiting systems. This would require a prohibitive number of Doppler-selected targets. In the 
absence of Doppler selection, a completely blind survey requires a more daunting average of 90 
nights per detection. 
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laDie i. 


Observing 


r Condit 


ions 




UTC date 




Seeing" 


Extinction'' ( 


mag) 


Mean PWV^ 






(") 


Median 


std dev 


95% value 


(mm) 


2010 Jun 22 


0.75 




0.187 


0.026 


0.238 


1.05 


2010 Jun 24 


0.78, 


1.07 


0.198 


0.045 


0.240 


1.13 


2010 Jun 25 


0.61 




0.180 


0.082 


0.241 


1.15 


2010 Jun 27 


0.99, 


0.74, 0.81, 0.8^ 


I 0.225 


0.197 


0.942 


2.48 


2010 Jun 29 


0.87, 


0.71 


0.237 


0.198 


0.894 


1.83 


2010 Jun 30 


0.68, 


0.73 


0.237 


0.202 


0.949 


1.43 


2010 Jul 27 


0.73 




0.213 


0.073 


0.319 


1.55 


2010 Jul 30 


1.02, 


1.03 


0.207 


0.093 


0.248 


0.85 


2010 Sep 15 


0.96, 


1.12 


0.215 


0.173 


0.922 


1.07 


2011 Jan 15 


0.83, 


0.88 


0.183 


0.056 


0.221 


1.74 



''From SNIPS focus images. 
^CFHT SkyProbe (V band) values. 

'^Based on mean of Caltech Submillimeter Observatory 225 GHz tau measurements. 



Table 2. Photometry Experiments 



UTC Date 


Experiments^ 


Integrations 


Cadence 


Aper. rad. 


Signal" 




rms** 






(sec) 


(min) 


(pixels) 


Target Reference 


xlO-4 


2010 Jul 27 


Snapshot (8-star) 


9 X 19 


17 


80 


15.8 




5.4 (5.8) 


2010 Jul 27 


Snapshot (8-star) 


9 X 31 


17 


75 


11.7 




6.0 (7.0) 


2010 Jul 27 


Snapshot (8-star) 


9 X 22 


17 


90 


15.7 




8.9 (5.2) 


2010 Jul 27 


Snapshot (8-star) 


9 X 30 


17 


75 


14.2 




7.3 (8.0) 


2010 Jul 27 


Snapshot (8-star) 


9 X 13 


17 


80 


11.8 




8.3 (9.0) 


2010 Jul 27 


Snapshot (8-star) 


9 X 23 


17 


85 


16.3 




7.8 (6.9) 


2010 Jul 27 


Snapshot (8-star) 


9 X 13 


17 


75 


13.6 




6.5 (7.5) 


2010 Jul 27 


Snapshot (8-star) 


9 X 18 


17 


75 


14.0 




8.4 (8.0) 


2010 Jul 31 


Common field 


15 X 88 


2.4 


65, 80 


8.4 


13.6 


9.7 (-) 


2010 Jul 27 


Common field 


13 X 22 


5.5 


60, 85 


13.1 


12.5 


8.4 (-) 


2011 Jan 15 


r Filter test 


20 X 50 


9.6 


85 


17.1 


16.8 


9.2 (9.5) 


2011 Jan 15 


z Filter test 


20 X 50 


9.6 


80 


18.0 


17.6 


8.5 (7.5) 


2011 Jan 15 


Exposure time tests 


22 X 25 


18.1 


65 


8.0 


7.1 


11.5 (11.7) 


2011 Jan 15 


Exposure time tests 


22 X 50 


18.1 


75 


16.1 


14.0 


9.8 (10.1) 


2011 Jan 15 


Exposure time tests 


22 X 100 


18.1 


90 


31.9 


27.8 


8.0 (11.5) 



^10^ e~ per integration. 

''Value in parentheses is theoretical precision (see §[5] and Table [3t . 

'^With the exception of the January 15 filter tests, all data listed here were taken with the z filter. 
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Table 3. Estimated Error/Noise Budget 



Noise Source" 


Median 


Mean 




10-4 


10-4 


Poisson 


4.8 


4.9 


Scintillation 


3.4 


3.3 


Atmospheric transparency fluctuations 


2.8 


2.9 


First-order extinction 


1.0 


0.91 


Sky 


0.99 


1.3 


Read and electronic pickup 


0.86 


0.93 


Motion on chip 


0.85 


1.0 


Second-order extinction 


0.45 


0.56 


Shape changes 


0.27 


0.34 


Chip/ Amplifier nonlinearities 


0.020 


0.021 



^Noise for each star is based on our model applied to each 
observation (theoretical and actual precision shown in Fig. I12I I. 
Noise is calculated per light curve (accounting for target and all 
comparison stars). 



-27- 




Earth around M4 star 



-Kepler (V=11.5) 
iCoRoT (V=11.5) 
J ground-based 



Neptune around MO star 




Jupiter around Sun 



-2-10 1 2 

hours after transit center 

Fig. 1. — Theoretical normalized light curves of stars transited by planets on 4 day orbits , comp ared 
with the photometric precisio n obtained from space - Kepler 30 minute (I Jenkins et al.ll2010l ) and 
COROT (jAigrain et al.l 12009) - and a p ractical limit from the ground (5 x 10^^). Limb darkening 
is included as described in lClaretl (2000|). A stellar mass-radius relation ~M;^ ''^ (appropriate for 
M dwarfs) is used. Ground-based photometry can be sufficiently stable to detect and characterize 
super-Earth- to Neptune sizes planets around M dwarf stars. 
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Hour Angle 

Fig. 2. — AitofF projection of trajectory of telescope pointings for constellation observations, de- 
signed to minimize slew time. The circle is approximately the size of the full Moon. Although a 

given target is usTially several degrees away from the next target, the slew time is not significantly 
longer than the readout time (^■20 s). The constellation technique enabled us to concurrently 
search for transits around eight target stars. 
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Wavelength (nm) 

Fig. 3. — Approximate transmission (bottom) and emission (top) of the atmosphere on a clear night 
over Mauna Kea with B, SDSS g, SDSS r, SDSS i, SDSS z, and narrow z filters multiplied by the 
SNIPS CCD quantum efficiency. There is a discontinuity at ~ 5200 A in the emission caused by 
low QE near the edge of the blue and red channels of the SNIFS integral field unit. Our observations 
were taken predominantly in SDSS z, although some were taken in SDSS r. Scattering affects the 
filters blueward of 5500 A whereas molecular (chiefly H2O and O2) absorption lines contaminate 
the r and z passbands. A narrow z filter (J. Johnson, private communication) can mitigate errors 
from both molecular absorption and scattering. 
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Fig. 4. — Photometric error in observations with the SDSS r (top left), SDSS z (top right), B 
(bottom left), and narrow z filters produced by a difference in spectral type between comparison 
and target star as a function of air mass (both using the SNIFS detector). The target star has an 
effective temperature of 4000 K and log (7 = 4.5. 
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Fig. 5. — Fractional error as a function of rms image motion within a bad region (high pixel 
response variation), a good region (low pixel response variation), and at the center of the SNIFS 
science-grade imaging detector. These curves were generated by taking a defocused image and 
shifting it around a map of the detector response produced from a series of ~225 fiats taken with 
varying exposure times through the SDSS r (top) and SDSS z (bottom) filters. Nonfilled points 
are for a pixel response map with a flat-field correction (median of 10 twilight flats). We assume 
no changes within the defocused PSF except position on the detector. Note that the good/bad 
regions used are comprised of different sets of pixels for the two different filters used. 
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Fig. 6. — Top: Mean of eight defocused images of a V = 13 star, each with an integration time of 
30 s. Bottom: Standard deviation of pixel values from the mean. All images have the same spatial 
scale (~ 8" on a side), but not the same intensity scale. The highes t peak in t he PSF is due to 
coma in the UH 2.2 m optics (C. Aspin, private communication 2010; lBehiill973l ). 



-33- 




100 200 300 400 

2 

X 

Fig. 7. — Histogram of the values of a fit to pairs of individual pixel values vs. median detector 
values in a series of 150 dome fiats. The theoretical distribution derived from a Monte Carlo analysis 
(based on Poisson statistics) is shown as a dashed line. We only include data from the bottom half 
of the SNIFS imaging science-grade detector (the region we use for transit imaging). Wc consider 
pixels with > 240 as poorly behaved and worth avoiding. These represent < 0.1% of the total 
number of pixels. 
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Fig. 8. — Comparison of photometric precision with and without fiat-fielding. Dome fiats are used 
when there is an insufficient quantity of high-S/N twilight flats from a given night. Although in 
some cases the precision is better with flats, on average, the precision is better without flat-fielding. 
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Fig. 9. — Curve of growth (vs. radius in arcseconds) of the mean of eight defocused images of a 
star obtained with the SNIFS/UH 2.2 m (sohd hne), and the S/N in an aperture of a given radius 
relative to the maximum value (dotted line). 
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Fig. 10. — Total predicted noise as a function of exposure time {top left) for a niz = 9 star. Also 
shown is the total noise {top right), optimal defocus radius {bottom right), and exposure time {bottom 
left), corresponding to the setup with the lowest theoretical noise as a function of z magnitude. 
We perform these calculations based on noise from Poisson statistics, read and electronic pickup 
noise, scintillation, sky background, and atmospheric transparency fluctuations. We have assumed 
1.5 minutes for readout and slew time between stars and that the comparison and target star are 
of identical flux and spectral type. Note that for our stars (mostly late K and early M dwarfs) 
V -z ~1.8 (jCovev et al. l boOTi '). 



-37- 



1.00 


1 1 1 1 1 1 1 1 1 1 
t » -•- 


1 1 1 1 1 1 1 
-•- * 


1 1 1 1 1 1 1 1 1 1 1 1 
~*~ ~ff~=~0.54x10"'~ 


1 1 1 1 1 1 1 1 1 
, -.- ^ 


1 1 1 1 1 1 1 1 1 1 1 1 
• - 
«- -• 

a 


1 1 1 1 1 1 1 
'=~1.07x10"'~~ 




• - r -•- ^ - • 


-•- • 


-•- - 


= 0.60x10^ 






= 0.72x10"^ 


0.98 




- , -•- ^ 


_ • _ 


~=~0.89x10^" 


• • 

»-«■- 




=1.30x10"^"" 








— • - 


~=~0.73x10"^~ 


-,- -9-^ -*- 




'=~0.83x10"^~ 


0.96 






• 


= 0.83x10"^ 


-•- ^- 




= 0.92x10"^ 




_^ ^_ ^_ 


» _•_ ^_ 




~=~0.78x10'^~ 


• • -•- 




'=~1 .72x10"^" 


0.94 








~=~0.65x10"^" 




• 


'="1.29x10"^"" 




" , , , 1 , , , 1 , , 


. 1 . . . 1 . 


-•- • 

. . 1 . 


~=~0.84x10"^~ 

, , 1 , , , 1 , " 


- ^ _•_ * _ • 

" , , , 1 , , , 1 , 


0- - 

• 

. 1 ... 1 ... 1 . . 


'=~^ .42x10^" 

. 1 . . . 1 . " 


0.00 0.02 0.04 


0.06 0.08 


0.10 


0.12 0.00 0.02 0.04 


0.06 0.08 0.10 


0.12 



JD - 2455404.769 JD - 2455404.769 

Fig. 11. — Light curves of eight stars using the constellation method along with rms using a 
reference signal constructed from stars observed immediately before and after the target (left), and 
from all seven other target stars (right). Cadence for these observations is 1 observation every 
17 minutes, which is sufficient to sample a transit (transits typically last > 1 hr). All light curves 
were normalized to 1, and then each was offset in increments of 0.01 for plotting. Although using 
more comparison stars improves the S/N, we achieve the best overall photometric precision for each 
target star using a reference signal built from the stars observed immediately after and before it. 



-38- 



1.2 



g 1.0 



0.8 



• • 



0.6 



_l I L 



0.6 



0.8 



1.0 



1.2 



Theoretical Photometric Stability (10" ) 



Fig. 12. — Precision of each of our observations vs. our estimated theoretical precision calculated 
for each observation. The dashed line marks where observed precision matches theoretical precision. 



4 6 8 10 

Doppler Threshold [m/s] 

Fig. 13. — Average number of photometric observing nights before a transit is detected in a snapshot 
survey of M dwarf stars with a Doppler signal above the specified threshold. In making this 
calculation, we take into account setup, slew, twilight, and integration time, and we use precision 
levels attained in our SNIFS observations. A planet population similar to those around G-type 
stars is assumed (see § 16.21 for details) 



